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Abstract. We present a protocol for increasing the entanglement between two 
entangled atomic ensembles based on applying an approximate atom- light beamsplitter 
transformation to both ensembles. The effective asymmetric atom- light beamsplitter is 
created via a double-pass quantum non-demolition interaction between polarized light 
and a spin polarized atomic ensemble, derived from the linearised dipole interaction. 
The entanglement concentration protocol itself uses the procrustean method, similar to 
that first devised for light by Browne et al [2003 Phys. Rev. A 67 062320] and includes 
photon counting after the interaction as the required non-Gaussian element. We 
calculate the output logarithmic negativity in this scheme and show that entanglement 
between macroscopic ensembles can be increased with probabilities comparable with 
those for the light scheme even if on-off detectors of low efficiency are used. 
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1. Introduction 

Quantum communication schemes, from quantum key distribution to teleportation and 
dense coding, offer better ways to exchange information in a network. Experiments 
into these areas are abundant and have yielded lots of successes over small distances. 
However, over larger distances, optical losses, phase diffusion, and mixing with thermal 
states cause the signals to decohere over some finite transmission length. In fact, the 
error probability scales exponentially with the length of the channel. 

The most obvious way to overcome this would be to create a series of signal 
amplifiers to be spaced out between the source and destination where the signal could 
be stored onto a quantum memory device and read out again. Such quantum repeater 
protocols have been devised yy with slight variances between schemes, but are in general 
reliant on entanglement distillation procedures and quantum memory devices that can 
be used to store a signal for a brief period of time. 

Due to the long lifetimes of atomic states, the most promising protocol so far uses 
the coherent spin states of cesium atoms to serve as a memory device. The quantum 
signal to be transmitted is contained in the polarization state of the incoming light mode 
and this information is written onto the macroscopic coherent spin state of the atomic 
ensemble. Julsgaard et al [2] successfully showed that the coherent spin states of two 
atomic ensembles could become entangled in such a way as to be analogous to the two 
mode squeezed state for light, written in the number basis 



where n is the photon number and A quantifies the reduction (squeezing) of the quantum 
uncertainty of the global state. We will discuss the particular meaning of these quantities 
in the case of atomic ensembles in Section [31 The better entangled the ensembles, the 
better the outcome of the quantum repeater protocol. However the entanglement of two 
atomic ensembles is at present limited and yields teleportation fidelities of < 0.6 [2J. For 
this reason, the option to increase entanglement in the ensembles by local operations 
would be greatly desired. 

Entanglement distillation [3l H] is a large field of research, and many different 
protocols have been devised for continuous variable entanglement concentration O |6l [7] . 
In recent years there has also been experimental success in distilling light states [HI |9l [10] 
but little is known about the possibility of distilling entanglement in atomic ensembles. 

We here demonstrate a theoretical protocol for increasing the entanglement between 
two atomic ensembles. Each atomic ensemble is made to interact with an incoming 
polarized light mode. Both emerging light modes are then redirected back into their 
respective ensembles and reemerge to be measured for reoriented photons using a photon 
counting technique implemented by on-off detectors. If both detectors respond in the 
affirmative, the entanglement between the two atomic ensembles can be shown to have 
increased. This scheme is similar in nature to one already devised for light (described 
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in [m [12] and the first demonstration reported in [8j). Note that the experimental 
schemes in [6l [71 [21 HO] present entanglement purification. That is, they begin with a 
mixed entangled state containing some non-Gaussian noise and aim at reducing that 
noise thus increasing the degree of entanglement. In contrast, the original procrustean 
scheme [HI HH [12] , the schemes based on the Kerr nonlinearity [5l [13] and the atomic 
scheme presented here deal with entanglement concentration. That is, the aim is the 
increase of entanglement content in an initially pure Gaussian partially entangled state. 

In Section [2]the quadrature system and interaction used in the scheme are described. 
In Section [3] the entanglement concentration is shown with losses taken into account in 
Section [H In Section [5] we conclude with discussion of the performance of the scheme. 



2. The system and interactions 

The standard QND Hamiltonian couples one of the two effective quadrature field 
operators of the strongly polarized light mode (re-scaled polarization variables) with 
one of the effective quadratures of the spin-polarized atomic ensemble, that is with a 
certain component of the re-scaled collective spin operator (for review see [13] )• The 
polarization variables are defined as follows. We assume the light beam to be strongly 
polarized in the x-direction so that the actual polarization variable, the Stokes operator 
Sx, can be replaced by its expectation value {S^). The Stokes operators can then be 
defined, e.g. for a pulse travelling in the z-direction, by 



c 



(alflx - alay) dr ^ (S^) = (2) 



T 





Sz = — [didy - aja^) dr (4) 
Jo 

where d^^y = d^^y {z, t) are the annihilation operators for photons linearly polarized in 
the X- and y- directions respectively with di {z, t) , a] {z', t) = 6ij6 {z — z') 6 (f — t') / c 

and is the real expectation value of d^ and d^ when highly polarized. T is the 
interaction time and t = t — z/c. With this in mind we can define continuous variable 
quadratures for the light state Xl and Pl as 



Here 



Xr.Pr 



-- i as with conjugate position and momentum and we have set h = 1. 
Note that these quadratures describe the quantum polarization of light and are not the 
usual amplitude and phase quadratures. 

A similar description for atoms is possible, also using appropriate, re-scaled 
"quadratures". Let us first define the relevant variables. The collective angular 
momentum of the atomic ensembles J, which has components Jj,j = x,y,z obey the 
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same SU(2) algebra as the Stokes operators, which makes it particularly easy to map 
quantum states of both systems onto each other. We are interested in increasing the 
entanglement of two macroscopic atomic ensembles, for example of two macroscopic 
ensembles of A^^^ cesium atoms at room temperature with a ground state degeneracy as 
used in Julsgaard et al [2]. There, in a homogenous magnetic field the cesium atoms 
were pumped into the \F = 4,mp = 4) state in the first cell and \F = A^nip = —4) 
in the second cell to form coherent spin states oriented in the +x and —x directions 
respectively. In this way, the collective angular momentum in the x-direction can also be 
replaced by the expectation value and (^Jxi^ = — (^Jx2^ with subscripts 1, 2 representing 
ensembles 1 and 2. In the language of density operators a^^^, = the operators can 

be written as ^ 

mp 

(7) 

mp 

Jz = — ^ (-^' "^-f) i^mp + l,mp — (^mp,mp+l) (§) 



2i 

mp 



where 



C (F, mp) = (F + l)-m^ (m^ + 1). (9) 



As the atomic ensemble is spin-polarized, conjugate position and momentum 
quadratures for atoms (subscript A) can be defined as (e.g.p^.[T5]): 

= Pa = -4^. (10) 

The interaction between light and atoms is represented by the linearised dipole 
interaction with far-off detuning (off- resonant interaction): 



^ = ^-dj-E(Rj) (11) 

i 

where dj = — erj is the dipole operator for the jth atom and Rj is the location of the 
jth atom. If, for example, the polarized light propagates in the z-direction through an 
atomic ensemble, the linearised interaction Hamiltonian can be written as 

H = a I S,{t)J,{t)dt ^ kPlPa (12) 
Jo 

where a is a coupling constant and k = {Sx) T {Jx)- Any higher order coupling terms 
are negligible if the laser beam is far detuned from the transition frequencies. 

Equation f[T^ is an example of a Quantum Non-Demolition (QND) Hamiltonian 
which would result in a phase shift in the Xl quadratures of the light by +kPa as it 
passes through the atoms. There is a corresponding back action on the distribution of 
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spins in the atomic ensemble, represented by a phase shift of +k,Pl in the Xa quadrature 
of the atomic ensemble. Physically, the interaction causes the polarization of the light 
to rotate about the axis of propagation, dependant on the quadrature distribution of 
the atoms. The back action effect on the atoms is to rotate the macroscopic spin state 
around the axis of propagation. 

Whereas a single pass of a light pulse through an atomic medium corresponds to 
the simple QND interaction described above, multiple passes open the possibility for a 
larger design freedom for the effective Hamiltonian, as for each pass a particular form of 
the underlying QND interaction can be adjusted (see e.g. [16]). A double pass scheme 
[17] can be used for the generation of polarization squeezed light by optical Faraday 
rotation. Another double pass scheme was suggested and thoroughly studied in the 
context of quantum memory for light modes based on macroscopic atomic ensembles at 
room temperature, as well as for the generation of entanglement between light and atoms 
[18j . There, the interaction Hamiltonian has been shown to include two main parts, 
one equivalent to the beamsplitter interaction, and the other to two-mode squeezing. 
Depending on the geometry of the setup, either of the two underlying dynamics can 
be selected. In this paper we exploit the fact that the double-pass Hamiltonian can 
approximate, with high fidelity, an actual beamspitter transformation, the quality of 
the approximation being dependent not only on the interaction strength, but also on 
the particular quantum states of the interacting modes [I9]. Tuning the interaction 
parameter, a highly asymmetric atom-light beamsplitter can be realized. In what 
follows, we describe the entanglement concentration scheme based on such an atom- 
light beamsplitter and analyse the performance of the scheme. 

3. Entanglement Concentration 

In the continuous variable regime, to increase the entanglement between two quantum 
objects with Gaussian quadrature distributions a non-Gaussian element is required. 
The measurement process offers this opportunity. Analogous to the photon subtraction 
schemes examined in [20] and ^ for increasing the entanglement in two mode squeezed 
states of light, a photon count heralds an increase in entanglement between two 
macroscopic atomic ensembles. 

In the case of the two mode squeezed state for light, the scheme is fairly simple, the 
two correlated beams are sent to separate highly transmissive beamsplitters, upon which 
they combine with another mode (usually the vacuum). There is a small probability 
that a photon will be subtracted from the main beam and proceed to be detected 
by a photon counter. If the counters at both beamsplitters register the presence of 
photons, then the entanglement of the two mode squeezed state increases. We wish for 
a similar "photon subtraction" scheme for our light and atoms in order to increase the 
entanglement between the two atomic ensembles as represented in Figure 1. 

The two atomic ensembles are prepared in an entangled two mode squeezed state 
as described in [2] . Previous works [21 [151 El] have detailed the methods for entangling 
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Figure 1. Procrustean entanglement concentration for atomic ensembles: the 
principle. The entanglement concentration protocol works similar to the light scheme 
[TT1[T2]. The two ensembles are entangled like a TMSS. Then a;-polarized light interacts 
with the atoms (number basis |0)) and non-Gaussian measurements are done to see 
whether the polarization of the photons has altered. A positive response from both 
detectors heralds an increase in entanglement between the atomic ensembles. The two 
square boxes symbolize the effective atom-light beamsplitter based on the double QND 
interaction. 




Photon 
Detector 



Figure 2. Procrustean entanglement concentration for atomic ensembles: details of 
atom-light interactions. Light is strongly polarized in the x-direction. It interacts 
with the spin-polarized atomic ensemble in a double-pass interaction equivalent to 
the effective atom-light beamsplitter which probabilistically accomplishes photon 
subtraction in a light beam. A y-polarized photon detection at the output implies 
that the atomic entanglement has increased. 



the collective atomic spins of the atomic ensembles using equations of the form ( |T2|) . By 
sending a polarized light mode through the atomic ensembles and taking a homodyne 
measurement, it is possible to collapse the atomic states into an entangled two mode 
squeezed state with variance A ^X^i — Xyi2 j = A (^Pai + pA-ij = where r is 
dependent on k and given by 

r = hn{l + 2K^). (13) 
Further, it is assumed that any displacement of the joint atomic quadrature distributions 
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from the centre in phase space caused by the entanglement process itself are small so 
that the constraints on the beamsplitter-like interaction still hold |19]. For ease of 
calculation, the atomic state is represented in the number basis as in Equation ([1]) 
where A = tanh(r) is the squeezing parameter, dependent on the interaction strength 
K between a light mode and atomic ensemble during the initial entanglement process 
via the relation f|T3|) . The number basis is considered to be the basis generated by the 
annihilation operator a a = + ipA up to normalisation. On consideration the number 
n represents the number of atoms in the upper excited spin state in the basis of the Jx 
operator. That is, n = 1 corresponds to the superposition of all possible combinations 
of atomic spins of the atoms in the ensemble for which a single atom is excited. 

If the light is initially prepared in a symmetric Gaussian state centred at = 0, 
Pl = then the light can be considered to be in a vacuum state. That is, there are no 
photons polarized in the 5"^ or Sz directions despite a steady base stream of linearly x- 
polarized photons so that (Sy) = {Sz) = 0. In the number basis generated by operators 
S+ = Sy + iSz = a/ {Sx){Xl + iPl) and S- = S\_ it is clear that a Fock state is simply 
the number of photons polarized linearly in the ^/-direction. 

In the number basis, the functionality of the photon subtraction scheme can be 
readily explained. If, after the beamsplitter-like interaction, one or more photons are 
detected behind the polarized filter, then the interaction has rotated the polarization of 
some photons in the light mode. The corresponding back action on the atomic ensembles 
is to flip the spin of one or more of the atoms in the system (see Fig. 2). 

The initial state of the light and atoms is given by 

oo 

= v/r^5^A"|n)i|n)2|0)3|0)4 (14) 

n=0 

where the first and second atomic modes are subscripted 1 and 2 respectively and 
the light vacuum modes are subscripted 3 and 4. The light vacuum modes 3 and 4 
(i.e. modes that are only polarized in the x-direction), are sent through the atomic 
ensembles 1 and 2 respectively and made to interact twice with the ensembles via the 
beamsplitter interaction described in [19]. That is, in one atom-light "beamsplitter" 
two QND interactions are performed. The first is of the form = (pP^XA and the 
outgoing light modes are redirected back into the atomic ensembles to interact via a 
second interaction, = —(pX^pA- We use the reasonable approximations on the 
quadratures that xi ~ a/1 — (pxi and (1 — 0^) Xa ~ a/1 — (pXa for small with similar 
conditions on pi and pa [19]. This double-pass scheme then performs the role of a 
beamsplitter transformation and is treated accordingly in what follows. The state after 
the beamsplitter interaction can be described by 

i..).^E^--i:/(55) 

n=0 fci,fc2=0 V \ V \ ^/ 

(1 _ ^2)2n-A,-fc2 ^ _ ^^^^j^ _ fc^) 2 1 ^1)3 1 ^2)4 (15) 

where is the strength of the interaction. 
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The outgoing light modes are then passed through a polarized filter to remove 
the base stream of x-polarized photons and allowing only photons whose polarizations 
have been rotated by the interaction until y-polarized to proceed (Fig. 2). A detector 
then registers whether y-polarized photons are present. As detectors are still relatively 
inefficient at counting photons, it is instead assumed that the detectors used can, to 
a high degree of efficiency, detect simply the presence of one or more photons ( "on-off 
detector"). The state thus becomes: 



f ^ 



S 

n=0 u,v=l 




X 0"+" (1 - ct^Y" " "\n- u)i\n - v)2 (16) 

(17) 

where S is the probability of getting an affirmative measurement at both detectors: 

1 - A2 2 (1 - A2) 



S = 



A2 (^2 + (1 _ ^2)2^ ' 1 - A2 (1 - 02)2 (^^2 + (1 _ ^2)2^ 



1 X2 

+ 4- (18) 

1-A2 (1-02)4 V J 

(19) 

The amount of entanglement in the two ensembles is quantified by the negativity and 
logarithmic negativity [22j of the state defined by 

m) = ^Tr ( v/(F^ - p^^) = , (20) 

£:^(p) = ln(l + 2A/-(p)) = ln(||p^^||) (21) 

where ||-|| denotes the trace-norm, i.e. \/ pp^, and p^'^ is the partial transpose of the 
density matrix p. 

It is simple to show that for a two mode squeezed state, the negativity and 
logarithmic negativity are given as 

Ar(TMSS) = (22) 
1 — A 

£:Ar(TMSS) = ln(l + A) -ln(l - A). (23) 

For projections onto an exact photon state (e.g., for a single photon detection), the 
negativity and logarithmic negativity can be calculated analytically. However, for on- 
off type measurements the entanglement measures must be calculated numerically due 
to the infinite sums over u and v. Fortunately, these sums appear to converge quickly 
and so we can set a reliable truncation point. The calculation is done in a similar way 
to the calculations of Kitagawa et al [2U]. Firstly, the density matrix of the state is 
expanded as 

oo 

|^A)out(^A| = Yl P«Ac,d|a)i(c| ® \b)2{d\ (24) 

a,b,c,d 



Entanglement concentration for two atomic ensembles using an effective atom-light beamsplitter9 
where 



Pa,b,c,d — 



i2\ °o 



u,v=l 



X 



a + [ a + u\ / c + u\ f c + u 
u / \ V / \ u / \ V 



X 



_ ^2)2(a+c+.-.) S^_^^^_J^_^^^_^_ (25) 

(26) 

The partial transpose of this state is given by 

oo 

{\'^A)out{^A\f^ = PaAc,b\a)i{c\ ® \b)2{d\ (27) 

a,b,c,d 

and the elements are zero unless the total Fock number of the entangled state, 
N = a-\-b = c + d, is non-zero. This follows from the delta functions in pa,d,c,b- 
Operator (I27|) is block diagonal and we can write it as a direct sum of each A^-dependent 
submatrix: 

il^AUi^Alf^ = ®^=o il^AUi^Alf^ (N) (28) 

where (|^A)out(*vi|)^^ (A^) is the (A^ + 1) x (A^ + 1) iV*^ submatrix. 
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Figure 3. A plot of the logarithmic negativity against A for (i) standard two 
mode squeezed state (solid) and for the output of the photon subtracted scheme 
for the beamsplitter like interaction between light and atomic ensemble when 
(ii) (p = 0.1 (dashed) and (iii) (p = 0.01 (dotted). See text for discussion. 
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Figure 4. Probability of success against A for (i) interaction strength (p = 0.1 
(solid) (ii) (j) = 0.05 (dashed) and (iii) (p = 0.01 (dotted). The probabilities of 
success are small but comparable with the light scheme. 

The negativity of the partially transposed state is then computed by numerically 
diagonalising each block individually to obtain the eigenvalues of each submatrix and 
adding up the absolute value of all negative eigenvalues. A cut-off, Nmax, must be 
introduced that is large enough compared to the mean number of excited spins. There 
is of course a trade-off between Nmax and the length of time needed to perform the 
calculation. For the purposes of this calculation, a value of N^ax = 100 gave very 
precise results. That is, at Nmax = 100 the numerical values of logarithmic negativity 
converge to at least 7sf and to increase Nmax beyond this is not beneficial. The results 
are shown in Figure 3. 

The entanglement between the two atomic ensembles is increased for all values 
of initial squeezing except for very high A. This increase in £j\f is not very large but 
comparable to when light modes are used in place of atomic ensembles [20], as are 
the probabilities shown in Figure 4. As can be seen, a trade-off is required between 
interaction strength and probability of success. As interaction strength increases, so 
does the probability of success, but the validity of the beamsplitter approximation 
decreases. Note, however, that the beamsplitter approximation has a very high fidelity 
of approximately 0.99 even for interaction strengths as high as ~ 0.35. For moderate 
initial squeezing, is largely unaffected by the interaction strength. As A approaches 
~ 0.95, the concentration procedure ceases to bring further benefit. However, this only 
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occurs at exceptionally high squeezing of the atoms, which is not experimentally viable. 
4. Modelling Detector Inefficiencies 

The largest contribution to loss in the photon subtraction scheme for hght modes comes 
from detector inefficiency. For the atomic ensemble scheme the efficiency of the detectors 
will also play a crucial role. The detectors here have a reduced number of photons 
to detect due to the polarization filter used to stop the base stream of x-polarized 
photons. The inefficiency of the detector can be modelled as an ideal detector behind 
a beamsplitter of transmittivity rj — p^. The light mode is combined with a vacuum 
on a beamsplitter and the vacuum is traced out before the projection measurement 
is performed. For a Fock state \k) combining with a vacuum, this amounts to the 
transformation 



1^' 0) ^ E J Q (vi^) ''^ \s, k 



s). (29) 

\/ \ .s / \ / 

Directly before the detection is performed, the state of the density matrix is given by 



,= (i-A^)f A- t ± ,/(:) (7) (;)(';) 

m,».0 fa,h=O«j,.0 V V*l/ \3lJ W \]2j 

X 0ii+fci+i2+fc2 _ ^2^2n+2m-ji-fei-j2-fc2 ^ ^ ^ ^ p^kiMJuji 

s=0 t=0 y=0 z=0 

\n - ki)i{m-ji\ (g)\n- A;2)2(m - J2I ® \s)3{t\ (g) \y)4{z\ (30) 



where 



A7-fel,fe2Jl,J2 

s,y,t,z 



X I/^ ^ ^dki-s,ji-t^k2-y,j2-z- (31) 

Light modes 3 and 4 are subsequently measured for the presence or absence of photons 
using the operator (1 — |0) (0|) and traced out. The density matrix of the two remaining 
atomic modes can then be described by 



min(m,n) 



CXJ 111111^, / , . 



+k2) 



m,n=Q ki,k2=0 
X (1 - ^2Yn+2m-2k,-2n Jl _ (1 _ _ (1 _ + (1 _ 

1^ — ki)i{m — ki\ ® |n — k2)2{fn — k2\ (32) 
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and the probability of success, taking into account detector losses, is given by 



5', 



1 - A2 [02 (1 _ ^) + (1 _ 02)2j + (1 _ 02)2j 

iT2- (33) 



1 - A2 [02 (1 _ ^) + (1 _ 02)2j 

(34) 

The effect that detector inefficiency has on the logarithmic negativity of the atomic 
ensembles can be shown in Figure 5. 



m — 





Figure 5. Dependence of the logarithmic negativity on the efficiency of the detectors. 
The soUd hne depicts the initial two mode squeezed state of the entangled atomic 
ensembles. Then the logarithmic negativity is shown for 77 = 1 (dashed line), 77 — 0.8 
(dotted line), 77 = 0.5 (dash-dot), and 77 — 0.2 (long dash). The interaction strength is 
= 0.1. 



As expected, the entanglement concentration becomes less pronounced for low 
detector efficiency. The positive message is that even for efficiencies as low as r/ = 0.2 
there is still a range of A values for which entanglement is increased (although the 
probability of success is quite low in this case). This range is experimentally accessible 
and so it is good news that the entanglement concentration protocol is more robust 
against imperfections. 
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5. Conclusions 

In this paper, we have employed a beamsplitter-like QND interaction between light and 
atomic ensembles based on the linearised dipole interaction between strongly polarized 
light and atomic levels for increasing the entanglement between two atomic ensembles. 
The entangled atomic ensembles in the initial two mode squeezed state interact with 
a highly polarized light mode in a quantum vacuum state via effective atom-light 
beamsplitter. The output light modes are subsequently detected using on-off detectors. 
This is analogous to the procrustean entanglement concentration scheme [HI |12] based 
on photon subtraction for distilling entanglement in light modes. Similar to their scheme, 
in our protocol detector clicks at both light outputs herald the successful "photon 
subtraction" (spin-flip in the atomic ensembles) and thus successful entanglement 
distillation. To assess the performance of the atomic entanglement distillation scheme, 
we have calculated the logarithmic negativity for the output quantum state of the 
two atomic ensembles and shown that it can increase indicating that the entanglement 
concentration procedure has been successful. The probability of success is very small 
but comparable to the probabilities for the corresponding light schemes that have been 
demonstrated experimentally [8]. The resulting atomic states are non-Gaussian and it 
remains to be seen whether they can be used for any teleportation procedures without 
devising a way to re-Gaussify the system. Use of the beamsplitter approximation allows 
us to closely mimic the entanglement concentration scheme that exists for light and to 
get an idea of how the light-atom protocol performs in comparison. The next step would 
be to remove the approximation and to assess the scheme for stronger interactions but 
a different approach has to be applied in this case. 
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